Kelvin Wave Cascade and Decay of Superfluid Turbulence 
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Kelvin waves (kelvons) — the distortion waves on vortex lines — play a key part in the relaxation 
of superfluid turbulence at low temperatures. We present a weak-turbulence theory of kelvons. We 
show that non-trivial kinetics arises only beyond the local-induction approximation and is governed 
by three-kelvon collisions; corresponding kinetic equation is derived. On the basis of the kinetic 
equation, we prove the existence of Kolmogorov cascade and flnd its spectrum. The qualitative 
analysis is corroborated by numeric study of the kinetic equation. The application of the results to 
the theory of superfluid turbulence is discussed. 

PACS numbers: 67.40.Vs, 03.75.Lm, 47.32. Cc 



The distortion waves on a vortex filament — Kelvin 
waves (KW) — have been known for more than a cen- 
tury P]. Superfluids with their topological (quantized) 
vorticity form a natural domain for KW 0. Nowadays 
there is a strong interest to the non-linear aspects of KW 
associated with studyinglow-temperature superfluid tur- 
bulence of *He H H H B i, H [i3L as well as vortex 
dynamics in ultra-cold atomic gase s |lltlT^. 

The superfluid turbulence 0,^3 is a chaotic tangle of 
vortex lines. In the absence of the normal component 
(T — > limit), KW play a crucial part in the vortex 
tangle relaxational dynamics. In contrast to a normal 
fluid, the quantization of the velocity circulation in a 
superfluid makes it impossible for the vortex line to relax 
by gradually slowing down. The only allowed way of 
relaxation is reducing the total line length. At T = even 
this generic scenario becomes non-trivial, as the total line 
length is, to a very good approximation, a constant of 
motion. In the scenario proposed by one of us Q, the 
vortex line length — in the form of KW generated in the 
process of vortex line reconnections — cascades from the 
main length scale (typical interline separation, i?o) to 
essentially lower length scales; ultimately decaying into 
phonons, as it was pointed out by Vinen 

A very specific feature of KW cascade is that the in- 
trinsic vortex line dynamics in the local-induction ap- 
proximation (LIA) (for an introduction, see, e.g., 
controlled by the small parameter 1 / ln(i?o / ^) , with ^ 
the vortex core radius, is subject to a specific curvature- 
conservation constraint rendering it unable to support 
the cascade process Q (see also below). Within LIA, an 
"external" ingredient of the vortex line dynamics — the 
vortex line crossings with subsequent reconnections — is 
required to push the KW cascade down towards arbi- 
trarily small wavelengths. The most characteristic fea- 
ture of this LIA scenario, distinguishing it from typical 
non-linear cascades, is the fragmentation of the vortex 
lines due to local self-crossings 3j; we will thus refer to 
this scenario as fragmentational scenario. 

Experimentally, the main consequence of the existence 
of a cascade regime, no matter what is its microscopic na- 
ture, is independence of the relaxation time of superfluid 
turbulence on temperature in the T — > limit. Davis 



et al. have observed such a regime set in in '^He at 
T < 70mK The numeric simulation of the vortex 
tangle decay at T = performed by Tsubota et al. 0| 
within the framework of LIA has clearly revealed the cas- 
cade regime; this may be considered as a circumstantial 
[l^ evidence of the proposed in Q scenario. 

A general question arises, however, of how far, in the 
wavenumber space, the structure of KW cascade is pre- 
determined by LIA dynamics. At wavelengths X Rq 
the non-local effects of the vortex line dynamics compete 
with LIA dynamics and can ultimately become the main 
driving force of the cascade. Moreover, given the specific 
spectrum of KW turbulence associated with the fragmen- 
tational scenario, where the amplitude of the turbulence 
is smaller than the wavelength only by a logarithmic fac- 
tor j^j, one concludes, that if non-local effects can in prin- 
ciple support the cascade, no matter how small is the cor- 
responding contribution at the main wavelength scale Rq, 
there will inevitably be such a wavelength scale A* <C -Rqj 
where the fragmentational scenario will be replaced by 
a purely non-linear — to be referred as pure — scenario in 
which vortex line self-crossings play no role. The exis- 
tence of the crossover between the two cascade scenar- 
ios has at least two important implications. First, the 
spectrum of sizes of vortex rings generated by decaying 
superfluid turbulence will have a lower cutoff A*. Sec- 
ondly, the spectrum of KW turbulence will be changed, 
which, in particular, is crucial for the cascade cutoff the- 
ory [H where the characteristic wavelength Aph at which 
KW essentially decay into phonons is a function of the 
cascade spectrum. 

A strong numeric evidence in favor of existence of pure 
KW cascade has been reported by Kivotides et al. 0]. 
Being very expensive numerically, this simulation was not 
able to accurately resolve the spectrum of KW turbu- 
lence. The data of recent simulation by Vinen et al. 
seem to be more conclusive. 

In this Letter, we propose a self-consistent (asymptot- 
ically exact in the limit of high wavenumbers) treatment 
for the pure KW cascade. More generally, we develop 
the KW kinetic theory in the regime of weak turbulence, 
where the smallness of non-linearities reduces the non- 
linear effects to scattering processes for the harmonic 
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modes. We find that the leading elementary process 
responsible for the kinetics is the three-kelvon scatter- 
ing. We demonstrate that the kinetics is entirely due 
to the spatially non-local interactions, the local contri- 
butions exactly cancelling each other. This cancellation 
sets the limit on the maximum power of the pure cas- 
cade: At the wavelength scale ~ Rq the contribution 
of the kelvon-scattering processes to the KW cascade 
contains a small factor l/ln^(_Ro/0 ^ compared to the 
reconnection-induced part. 

In terms of weak turbulence theory, KW cascade is a 
Kolmogorov cascade associated with the transport 
of energy (closely related to the vortex line length in our 
case) in the wavenumber space. We establish the Kol- 
mogorov spectrum and find the relation between the en- 
ergy flux and the amplitude of KW turbulence. Finally, 
on the basis of our results we discuss the fragmentational- 
to-pure cascade crossover, which we predict to be rather 
extended (up to two decades) in the wavenumber space. 

Hamiltonian. We employ the Hamiltonian representa- 
tion of the vortex line motion Q, which is exact up to 
a certain geometric constraint: there should exist some 
axis z, with respect to which the position of the line can 
be specified in the parametric form x = x{z), y = y{z), 
where x and y are single- valued functions of the coordi- 
nate z. This perfectly suits our purposes, since we are 
interested in the wavelengths substantially smaller than 
i?o a-nd thus will treat the vortex line as a straight line 
with small-amplitude distortions. In terms of the com- 
plex canonic variable w{z, t) = x{z, t) -\-iy{z, t), the Biot- 
Savart dynamic equation acquires the Hamiltonian form 
iw = SH[w]/5w*, with ^ 



K f, , 1 + Rew'*{zi)w'{z2) 

47r7 ^(^zi- Z2Y + \w{zi)-w{z2)\'^ 



(1) 



where k = 2Trh/m is the circulation quantum (m is the 
particle mass). Our fundamental requirement that the 
amplitude of KW turbulence be small as compared to 
the wavelength is formulated as 



\w{zi) - W{Z2)\ 

a[zi,Z2) = ; : < 1 



\Zl - Z2\ 



(2) 



This allows us to expand in powers of a: H ^ + 
Hq + Hi + H2 + ■ ■ ■ {Eq is just a number and will be 
ignored). The terms that prove relevant are: 



Hn = — 



dzidz2 

\Z1 - Z2\ 



[2Rew'*{zi)w'{z2)~ a^] , (3) 



Hi 



Ho = 



327r7 \zi-Z2\ 



647r J \zi — Z2I 



The Hamiltonian Hq describes the linear properties 
of KW. It is diagonalized by the Fourier transformation 



w{z) = L flfcc*'^^ {L is the system size, periodic 

boundary conditions are assumed): 

Ho = ^J2''k afeflfc , i^k = («:/47r) Hl/k^) , (6) 

A; 

yielding Kelvin's dispersion law ujk ■ 

Though the problem of KW cascade generated by de- 
caying superfluid turbulence is purely classical, it is con- 
venient to approach it quantum mechanically — by intro- 
ducing KW quanta, kelvons. In accordance with the 
canonical quantization procedure, we understand ak as 
the annihilation operator of the kelvon with momentum 
k and correspondingly treat w{z) as a quantum field. A 
minor caveat is in order here. The Hamiltonian func- 
tional 1^ is proportional to the energy — with the coeffi- 
cient Kp/2, where p is the mass density — but not equal 
to it. This means that if one prefers to work with gen- 
uine Quantum Mechanics rather than a fake one (for our 
purposes, the latter is also enough), using true kelvon an- 
nihilation operators, a^j, field operator w, and Hamilto- 
nian, H, he should take into account proper dimensional 
coefficients: H = {kp/2)H, w{z) = \/2h/KpLJ2k Ofce''^^. 
By choosing the units ?i = K = l,p = 2, we ignore these 
coefficients until the final answers are obtained. 

In the quantum approach, there naturally arises the 
notion of the number of kelvons. This number is con- 
served by the Hamiltonian in view of its global U{1) 
symmetry, w — > e'^^w (reflecting rotational symmetry of 
the problem). Another advantage of the quantum lan- 
guage in weak-turbulence problems is that the collision 
term of the kinetic equation immediately follows from the 
Golden Rule for corresponding elementary processes. 

Kelvon scattering. As was mentioned above, Hamil- 
tonian implies only elastic scattering. The one- 
dimensional character of the problem in combination 
with the conservation of the momentum and energy 
suppresses the two-kelvon scattering channel: the pro- 
cess (fci,fc2) {k^jki) is only possible if either {k^ = 
ki, /c4 — k2), or (fca — k2, k^ — ki) which does not 
lead to any kinetics. We thus conclude that the lead- 
ing process in our case is the three-kelvon scattering. 
The processes involving four and more kelvons a much 
weaker due to the non-equality (|2Jl. The effective ver- 
tex, Vi2f, for the three-kelvon scattering process [sub- 
scripts (superscripts) stand for the initial (final) mo- 
menta of the three kelvons; we use a short-hand nota- 
tion replacing each momentum kj with its index j] con- 
sists of two different parts. The first part involves terms 
generated by the two-kelvon vertex, A (corresponding to 
the Hamiltonian Hi ) in the second order of perturbation 
theory. [All these terms are similar to each other; we 
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explicitly specify just one of them: A{2G{oj7,k7) A^' ^. 
Here G{uj,k) = l/{uj — uJk) is the free-kelvon propaga- 
tor, uji = uji + LU2 — uj^, fcy = fci + fc2 ~ k^.] The second 
part of the vertex V is the bare three-kelvon vertex, B, 
associated with the Hamiltonian i?2- 

The explicit expressions for the bare vertices directly 
follow from Q and ©: A = {6D - E)/8tt, D'^f l = 
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/(f (da;/x5)[l-C{i}-C{2}-C^3}-CW+c|2j + C{43>+ occupation numbers terms: 



G 



{4}i 
{2}J 



Cl^'} ) + A;3fci ( + - C{43} _ cgj ) + fc3fc2( Ci3} + 
B = (3P- 5g)/4^, Pi%%^ = /o''(da;/x5) fc6fc2[C{2} - 

'-"{2}~*-'{23}+<-'{23}~'-^{2}+'-'{2} + ^-"{23} ~ '"'{1} ~ 
Ci56} _ Cgj + eg} _ ^{46} ^ ^{456} ^ c-f^e} _ ^^^^^j^ 



4,5,6 
2,3 



_ c<{456} _ ^{56} + C{23} ^ ^3} + Cg + ^{13} 



'{13} 
^{5} 



'{3} 
^{6} 



c 



{46} 



{3} 
(^{65} 



^{13} ^{2} 



(^{5} _ (-,{45} 



{3} ^{3} 

{13} ^{2} ~ ^{3} ^ "^{la} + <^{2} ~ <^{2}]- Here C's 
denote similar looking cosine functions: C^ij ~ cosfcix, 

fcs - fci)x, 



clt} — cos(A;4 



"{i} 

{45} 



ki)x, C|^|^ = cos(A;4 



{12} ~ cos(fc4 + fcs — /ci — fc2)2^7 and so forth. The in- 
tegrals for bare vertices A and B logarithmically diverge 
at their lower limits and, generally speaking, one has to 
introduce a cut-off parameter ^. However, in the fi- 
nal answer for the effective vertex V the logarithmic in 
^ terms cancel each other (it is important here to take 
into account that, in the propagator G, the frequency LUk 
also diverges logarithmically with ^ ^ 0), so that V is 
finite in the limit ^ — > 0. This fact is not just a mere 
coincidence. The small- a; contributions to the integrals 
correspond to the local- induction dynamics, which can- 
not lead to any kinetics because of the extra constants of 
motion Cancellation of the LIA contributions at the 
level of the effective scattering amplitude is an explicit 
demonstration of this circumstance. 

Kinetic equation. The kinetic equation is written in 
terms of averaged over the statistical ensemble kelvon 
occupation numbers Uk = {al,ak)- 



1 



ni 



(3-1)! 3! 



/4,5,6 
,2,3 



(7) 



Here W^2^'3 is the probability per unit time for the 
elementary scattering event (^1,^2,^3) — > (^4,^5,^6); 
combinatorial factor compensates multiple counting the 
same scattering event. For our interaction Hamilto- 
nian Hi^t = J2ki fee ^(^^) ^1*2^3^ a6a5«4i3Q'2'^i [where 
the vertex V is obtained from V by symmetriza- 
tion with respect to corresponding momenta permu- 
tations; S{k) is understood discretely as Sk,o, and 
Afc = ki + k2 + k^ — k4 — k^ — fcg] the Golden 
Rule reads: <2%^ = 2^|(3!)2y,^-5'6|2^4,5.6 ^^^^^^^^^^^ 

where /f'2.'3 = nin2nz{ni + l)(n5 -f l)(n6 -|- 1), Aw = 
wi -|- a;2 + W3 — a;4 — — cje- Combinatorial factor 
(3!)^ takes into account addition of equivalent ampli- 
tudes. The classical- field limit of the quantum kinetic 
equation Q is obtained by retaining only the largest in 



fii— 2167r 



fc2,...,fc( 
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(8) 

where 23 = 'T'i'^2f^3(f^4'T-5 + n^riQ + n^ng). 

Kolmogorov cascade. Kinetic equation ||SJ) supports 
Kolmogorov energy cascade jl5,|, provided two condi- 
tions are met: (i) the kinetic time is getting progressively 
smaller (vanishes) in the limit of large wavenumbers, (ii) 
the collision term is local in the wavenumber space — not 
to be confused with the local-induction approximation 
in the real space, — that is the relevant scattering events 
are only those where all the kelvon momenta are of the 
same order of magnitude. We make sure that both con- 
ditions are satisfied in our case: the condition (i) can be 
checked by a dimensional estimate, provided (ii) is true. 
The condition (ii) is verified numerically. 

The cascade spectrum can be established by dimen- 
sional analysis of the kinetic equation. The estimate of 
Eq. (IHl yields: iik ^ k^ ■ \V\^ ■ uj^^ ■ k~^ ■ nl, the factors 
go in the order of the appearance of corresponding terms 
in ((SJ. At ki . . . ke k we have \V\ ^ k^ and 



nk 



(9) 



The energy flux (per unit vortex line length), 9k, at the 
momentum scale k is defined as 



L ^ ^ cjfc/ ilk' 

k'<k 



(10) 



implying the estimate 9k ^ khkLOk- Combined with ||^, 



this yields 9k 



51,17 



nlk 



, and the cascade requirement that 



9k be actually /c-independent leads to the spectrum: 

(11) 



nk 



The value of the parameter A in (|ll|l corresponds to the 
value of the energy flux. An accurate relation between 9 
and A is (we restore all dimensional parameters): 



61 « 3 • 10" 



fi^A^ 

K2p4 



(12) 



The dimensionless coefficient in this formula has been 
established numerically (see below) with the error ~ 75%. 
One should not be confused with the appearance of the 
Plank's constant in this relation — it is entirely due to the 
fact that we use quantum mechanical definition of n^. 

The exponent 17/5 is in a perfect agreement (within 
the error bars) with the spectrum observed in Ref. 
{Q. of Ref. JJl] is proportional to our rife); it fits the data 
much better than the exponent 3 used by the authors. 

It is useful to express the KW cascade spectrum 
in terms of a geometric characteristic — typical ampli- 
tude, of the KW turbulence at the wavevector ^ 
k. By the definition of the field w{z) we have: 6^ ^ 
Eo^fc(4"9) = T.q^k ^ ^"fe- Hence, 



hk oc ^i/2A:-6/5 



(13) 



4 



One can convert Hll() into the curvature spectrum. For 
the curvature c{() — d^s/dC,^ [where s(C) is the radius- 
vector of the curve as a function of the arc length C,]^ 
the spectrum is defined as Fourier decomposition of the 
integral Ic — f |c(C)P dC- The smallness of a, Eq. (0), al- 
lows one to write Ic ~ J dz{'w"^z)w" {z)) — fc^ Uk oc 
fc'^/^, arriving thus at the exponent 3/5. 
Numerics. The aim of our numeric analysis is (i) to 
make sure that the collision term of the kinetic equation 
is local and (ii) to establish the value of the dimension- 
less coefficient in II12II . The analysis is based on the fol- 
lowing reasoning . Consider a power-law distribution 
of occupation numbers, Uk — A/k^, with the exponent 
(3 arbitrarily close, but not equal, to the cascade expo- 
nent Po = 17/5. Substitute this distribution in the colli- 
sion term of the kinetic equation — right-hand side of ||SJ). 
Given the scale invariance of the power-law distribution, 
the following alternative takes place. Case (1): collision 
integral converges for /3's close to /3o, and, in accordance 
with a straightforward dimensional analysis, is equal to 



Coll(K - A/k'^],k) = C{(3)A^uj-^^/k^''^~ 
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(14) 



Here C(/3) is a dimensionless function of /3, such that 
C(/3o) — since the cascade is a steady-state solu- 
tion. Case (2): collision integral diverges for /3 close 
to Pq. The case (2) means that the collision term is 
non-local and the whole analysis in terms of Kolmogorov 
cascade is irrelevant. [Fortunately, our numerics show 
that we are dealing with the case (1).] Substituting 
((njl for hk in ((Tn|) . we obtain the expression 9k{(3) = 
(y4V27r)fc"-5''C(/3)/(17-5/3). Taking the limit /3 ^ /3o, 
we arrive at the /c-independent flux 9 = —C'{(3o)A^/10tt. 
We use this formula to obtain the coefficient in p2|l by 
calculating C(/3) and finding its derivative C'{Po). We 
simulated the collision integral by Monte Carlo method. 
The integrals D, E, P,Q, were calculated numerically. 
The slowing down of the simulation due to the latter 
circumstance is the main source of large relative error. 



Superfluid turbulence. The main characteristic of the 
superfluid turbulence is the vortex line density, L (to- 
tal line length per unit volume; note that Rq ^ 1/VL.) 
The parameter A is related to L through the energy 
flux 9. Since (to a good accuracy) the energy is the to- 
tal line length times {pn"^ /An) ln(i?o/f), the same coeffi- 
cient relates the energy flux to the line-length flux. The 
line- length flux, L/ L, is available from the simulation of 
Ref. 0: L/L « 6 • IQ~^kL ln(l/eVI). We thus get 



{fiA/npf « 2L In2(l/^VI) 



For the amplitude spectrum this yields: 
hkk - [(\/I/fc) ln(l/CVT) 



1/5 



(15) 



(16) 



With this spectrum, Vinen's prediction for the cutoff 
momentum fcph ^ X^^^ based on Eq. (2.24) of Ref. 
should read (c is the sound velocity): 



L 



ln^/^(l/^^/Z) 
lni°/^(l/^A:ph) 



-1 5/6 



(17) 



From it is seen that the crossover from the frag- 
mentational cascade regime to the pure one should be 
very slow. Indeed, to significantly suppress the fragmen- 
tational regime — that is to suppress local self-crossings 
of the vortex line — it is necessary to make the parame- 
ter hkk substantially smaller than unity. In view of the 
exponent 1/5, this requires increasing fc by ~ two or- 
ders of magnitude. In this extended crossover region, one 
has hkk ~ 1, while for the applicability of the treatment 
developed in this Letter one needs hkk <^ 1. Hence, a 
reliable description of the crossover from the fragmenta- 
tional to the pure cascade seems to be impossible without 
a direct numeric simulation of vortex line dynamics. 

We are grateful to S. Nazarenko, W. Vinen, 
C. Barenghi, and M. Tsubota for useful discussions. 
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